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ABSTRACT 

The application of the virial theorem to the Broad Line Region of Active Galactic Nuclei allows Black 
Hole mass estimates for large samples of objects at all redshifts. In a recent paper we showed that ionizing 
radiation pressure onto BLR clouds affects virial BH mass estimates and we provided empirically calibrated 
corrections. More recently, a new test of the importance of radiation forces has been proposed: the Mbh-o" 
relation has been used to estimate Mbh for a sample of type-2 AGN and virial relations (with and without 
radiation pressure) for a sample of type-1 AGN extracted from the same parent population. The observed 
L/LEdd distribution based on virial BH masses is in good agreement with that based on Mbh-o" only if radiation 
pressure effects are negligible, otherwise significant discrepancies are observed. In this paper we investigate the 
effects of intrinsic dispersions associated to the virial relations providing Mbh, and we show that they explain 
the discrepancies between the observed L/LEdd distributions of type-1 and type-2 AGN. These discrepancies 
in the L/Lndd distributions are present regardless of the general importance of radiation forces, which must 
be negligible only for a small fraction of sources with large L/Lndd- Average radiation pressure corrections 
should then be applied in virial Mbh estimators until their dependence on observed source physical properties 
has been fully calibrated. Finally, the comparison between Mbh and L/LEdd distributions derived from -based 
and virial estimators can constrain the variance of BLR physical properties in AGN. 

Subject headings: radiation mechanisms: general — galaxies: active — galaxies: fundamental parameters — 
galaxies: nuclei — quasars: emission lines — galaxies: Seyfert 



1. INTRODUCTION 

In the last few years, it has become increasingly clear that 
supermassive black holes (BH) are an essential element in the 
evolution of galaxies. The key observational evidence of a 
link between a BH and its host galaxy is provided by the tight 
correlations between BH mass and luminosity, mass, velocity 
dispersion and surface brightness profile of the host spheroids 



(e.g. IGebhardt etal.'2000', 'Ferrares e & MerrittpOOO] [Marconi 
[S Hunt 2003 ; Graham & Driver 200]}. The link between BH 
and host galaxy is probably established by the feedback of 
the accreting BH on the host galaxy itself (e.g. [Granato et al. 
[2004l|Di Matteo et al.|2005l|Croton et al.|2006| and references 
therein). 

In order to fully understand the implications of BH growth 
on host galaxy evolution it is essential to measure BH masses 
in large samples of galaxies from zero to high redshifts. Since 
direct BH mass measurements based on stellar and gas kine- 
matics are possible only in the local universe (e.g. Ferrarese 
|& Ford] [2005 | l, less direct estimators have been calibrated 
following a "BH mass ladder" (Peterson 20()4'). The final 
rung of this ladder is provided by the virial estimators which 
allow us to estimate BH masses from the spectra of AGN 

Electronic address: marconi@ai'cetri. astro. it 

' Dipartimento di Astronomia e Scienza dello Spazio, Universita degli 
Studi di Firenze, Largo E. Fermi 2, 50125 Firenze, Italy 

^ Physics Department, Rochester Institute of Technology, 85 Lomb Memo- 
rial Drive, Rochester, New York 14623, USA 

' INAF - Osservatorio Astronomico di Roma, Via Frascati 33, 1-00040 
Monte Porzio Catone, Italy 

Research Center for Space and Cosmic Evolution, Ehime University, 2-5 
Bunkyo-cho, Matsuyama 790-8577, Japan 

' INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125, 
Firenze, Italy 

* Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cam- 
bridge MA 02138, USA 



with broad emission lines (type-1): under the assumption that 
the Broad Line Region (BLR) is gravitationally bound and 
its motions virialized. The BH mass can be expressed as 
^BH = fV^ Rblr/G, where Rblr is the BLR average distance 
from the BH, V is the width of the broad emission line and / is 
a scaling factor which depends on (unknown) BLR properties. 
i^RT R can be estimated with the Rms -L relation (Rblr oc L"^, 
Kaspi et al.| 2000','Bentz et al.'2008^ leading to Mbh = fV^L^ 
where / is calibrated empirically st arting from the Mbh ~ c re - 
lation (e.g. Onk en et al.|20 04; Vest ergaard & Peter son'2006'). 

One of the basic assumptions of reverberation mapping is 
that the BLR is photoionized (Blandford & McKee 1982) im- 
plying that BLR clouds are subject to radiation forces aris- 
ing from ionizing photon momentum deposition. In a recent 
paper (Marconi et al.1'2008 , hereafter M08), we showed that 
these radiation forces constitute an important physical effect 
which must be taken into account when computing virial BH 
masses. We empirically calibrated a radiation pressure correc- 
tion of the form Mbh = fy^L'' + gL, and we showed that it is 
consistent with a simple physical model in which BLR clouds 
are optically thick to ionizing radiation and have average col- 
umn densities of A^h — lO^"* cm"^ towards the ionizing source. 
This value is remarkably similar to that adopted in standard 
photoionization models to explain observed BLR spectra. 

Recently, ["Netzer| (|2009) (hereafter N09) proposed a test of 
the importance of radiation forces on virial BH mass esti- 
mates. He selected two large samples of type-2 and type-1 
radio quiet AGN drawn from the SDSS survey, and covering 
the same range of redshift (0.1 < z < 0.2) and continuum lu- 
minosity (10^^-^ < ALa(5100A) < 1044-8 ergs"'). After elimi- 
nating type-2 galaxies classified as LINERS resulted in a final 
sample composed of 4197 and 1331 in the type-2 and type-1 
objects respectively. By comparing the distributions of [O III] 
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line luminosities he concluded that the two samples were ex- 
tracted from the same parent population. After the sample 
selection, N09 estimated BH masses for the type-2 sample us- 
ing the Mbh - c relation in the calibration by Tremaine et al.| 
(|2002 ). For the type-1 sample, he estimated BH masses using 
both the classical virial relation (e.g. jVestergaard & Peterson 



2006 1 and the one by M09 with the correction for radiation 
pressure. He then compared the distribution of BH masses 
and L/Lem ratios finding significantly different distributions 
between type-2 and type-1 AGN under the assumption that 
radiation pressure significantly affects BLR motions. In par- 
ticular, while the L/Lndd distribution of type-2 AGN is broad, 
bell-shaped and extends up XoL/Lem ^ 1, the L/Leu distribu- 
tion of type-1 AGN with an important radiation pressure cor- 
rection is strongly peaked at L/ Lndd — 0.15 with a sharp cutoff 
at L/Lndd ^ 0. 15 (see fig. 3 of N09). Conversely, the distribu- 
tion of type-1 AGN without the radiation pressure correction 
is in good agreement with that of type-2 AGN. N09 then con- 
cluded that radiation forces are not affecting BLR motions, 
which is possible only if BLR clouds have extremely large 
column densities (A^h ^ 10^'*cm"^). 

In this paper we generalize the conclusions by N09 and we 
investigate the origin of the discrepancies amongst the L/L^m 
distributions. The observed discrepancies do not imply that 
radiation forces are not important. They are expected from 
the physically-justified mathematical expression of radiation 
pressure corrected virial masses, if the intrinsic dispersion as- 
sociated to such scaling relations is not taken into account. In 
particular, to obtain the "true" Mbh and L/Lem distributions 
one should take into account the dispersion of the scaling pa- 
rameters used in the virial estimators. We will not address the 
issue as to whether outflows affect BLR motions (e.g. Chiang] 
|& Murray! 1996[ [Kurosawa & Proga|2008/ but we will assume 
that BLR clouds are gravitationally bound. 

In ^we discuss why truncated L/Ledd distributions are ex- 
pectedwhen using virial mass indicators with the correction 
for radiation pressure. In ^we use Montecarlo simulations 
to reproduce the results of N09 and we show that similar ob- 
served distributions are expected regardless of the importance 
of radiation pressure. Finally, we summarize our results and 
draw our conclusions in ^ 

2. OBSERVED AND TRUE DISTRIBUTIONS OF EDDINGTON 
RATIOS 

In this section we provide the physical explanation of the 
expected differences between "observed" and "true" distribu- 
tions of Mbh and L/Lndd- For a given physical parameter, we 
use "true" to denote its distribution of true values which is not 
actually observable because of measurement errors or uncer- 
tainties in adopted scaling relations. We use "observed" to de- 
note the distribution of values obtained from the observations, 
by direct measurements or by applying scaling relations. 

The classical version of the virial theorem which does not 
take into account radiation pressure provides a mass estimator 
which can be written as ( [Vestergaard & Peterson|2006 i: 
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BH = 



lO^y^^L^^iooMo 



(1) 



where / is a scaling factor which encodes BLR geometry, 
physical structure and projection effects, Vnp is the FWHM 
(Full Width at Half Maximum) of the broad H/3 line (units of 
1000 km s"'), and L5100 is the continuum luminosity at 5100 A 
{XLx, in units of 10'*'* ergs"'). For consistency with N09 we 
adopt / = 6.7 and slope 7 = 0.6. When taking into account 
radiation forces due to the absorption of ionizing photons, the 



virial mass estimator is modified as follows (M08): 

Mbh = 10-^"'^ yJ^L^iooMo + lO^LjiooM© (2) 

The second term represents the correction of radiation pres- 
sure and, in the assumption that BLR clouds are optically 
thick to ionizing photons and neglecting other sources of ra- 
diation pressure, it is given by (MOB): 
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where Lion is the AGN luminosity in H-ionizing photons, A^h 
is the total column density of BLR clouds (of both ionized 
and neutral gas) towards the ionizing source, bion and b^im are 
the bolometric corrections for the ionizing and optical contin- 
uum luminosities, respectively (Lboi = ^ionAon = ^sioo^sioo)- 
The empirical calibration performed by M08, and adopted by 
N09, assumed 7 = 0.5 and provides /,-ad — 6.13 and g' ~ 7.72; 
the latter value indicates an average BLR column density 
A^h = 10^^ cm"^, consistent with expectations from photoion- 
ization models. The column density of the BLR clouds sets 
the relative importance of the gravitational force (which de- 
pends on cloud mass) and the radiative force (which is inde- 
pendent of mass). Therefore it is the most critical parameter 
for the radiation pressure correction and we can outline its 
effects by writing: 
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where A^23 is A^h in units of 10^^ cm ^. The Eddington ratios 
with and without radiation pressure correction are therefore: 
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where LEdd.o is the Eddington luminosity for a 1 Mq ob- 
ject. The observed distribution of Eddington ratios in the pure 
virial case is determined by the observed distributions ofVHp 
and Lsioo- For instance, if Vnp and Lsioo are log-normally 
distributed, then so are the Eddington ratios. Similarly, the 
observed distribution of Eddington ratios with the radiation 
pressure correction is determined by the observed distribu- 
tions of Vhi3 and L5100 but with an important additional fea- 
ture which results from the addition of the radiation pressure 
correction. When the radiation pressure correction dominates 
over the virial term, ^0^""'~^V^pL]^^N23 ^ 1, the Eddington 
ratio becomes asymptotically constant 
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where Lblr(A^h) is the critical luminosity at which radiation 
forces on BLR clouds balance gravitation, and is a function of 
Nh (Lblr A^h)- The physical meaning of this behaviour is 
straightforward: the virial mass estimator with the correction 
for the radiation pressure is valid only if the BLR is grav- 
itationally bound, that is if Lboi < Lsu^iNn) < L^dd- Since 
LuLRiNH) is always smaller than Lsdd, LBLRiNu) / LEdd is the 
maximum allowed Eddington ratio and any L/Lndd distribu- 
tion derived from eq.[4|will be truncated at that value. 

At first sight it might be assumed that the observed distribu- 
tions of Lboi/LEdd, obtained from eqs.|5]and|6]and depending 
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only on the observed distributions of Vhp and L^ioo, provide 
a good description of the "true" ones. However this propo- 
sition does not take into account the fact that the scaling pa- 
rameters themselves (/, /i-ad, g, A^23, ^5100) are drawn from 
their own true distribution functions and necessarily vary from 
one object to another. In order to properly represent the true 
Lboi/^Edd distribution this dispersion must be included, and 
the interpretation of the observed differences between true 
and observed Mbh and Lboi/^Edd distributions must be mod- 
ified accordingly. Since the scatter in the scaling parameters 
is not known and cannot be taken into account, the observed 
Lboi/i'Edd distribution in the case of the radiation pressure cor- 
rection will always present a sharp cutoff at iBLR/i'Edd- This 
cutoff is smeared away when taking into account the true dis- 
tribution of scaling parameters in general, and of BLR col- 
umn densities in particular. |Onken et al.| ( |2004| l showed that 
the r.m.s. scatter of the Mbh-ct relation using virial masses 
is ~ 0.5 dex, compared to the ~ 0.3 dex of the same relation 
using more direct BH mass measurements from spatially re- 
solved stellar or gas kinematics. Part of this additional scatter 
is probably explained by a broad distribution of / (or /i-ad) val- 
ues, which is naturally expected since the physical properties 
of BLR clouds must be characterized by a variance from one 
object to another (e.g., different cloud geometries and spatial 
distributions, relative orientations of the line of sights, etc.). 
Recently, fast eclipsing of the X-ray emitting source in the 
Seyfert galaxy NGC 1365 has been unambiguously explained 
by occultation from fast moving BLR clouds ( |Risaliti et al. 
|2007||2009] l. In particular, these fast eclipses allow us to esti- 
mate the column density of BLR clouds toward the AGN, A^h 
finding a distribution of values in the 10^''-10^'*cm"^ range. 
The existence of a relatively broad distribution of A^h values 
in a single object might indicate an even broader distribution 
of values over the whole AGN population. 

The L/^Edd distribution plotted in Figure 3 of N09 clearly 
shows the features discussed above. In particular, the dis- 
tribution obtained with the correction for radiation pressure 
is sharply cutoff at L/Lem — 0.15 which corresponds to the 
critical value at which the BLR becomes gravitationally un- 
bound. Conversely, the L/Lndd distributions for type-2 AGN 
and type-1 AGN with classical virial masses are bell-shaped, 
with a tail beyond L/L^dd — 0.15. N09 interpreted these 
differences in the L/^Edd distributions as an indication that 
BLR clouds have column densities significantly larger than 
lO^^'cm"^, for which the radiation-pressure force term is neg- 
ligible. However, a distribution of relatively small A^h values 
(~ lO^-'cm"^) with a tail extending to large values could in 
principle provide an alternative explanation of the existence 
of type-2 AGN with L/Leu beyond the critical value, without 
implying that radiation forces are not important in general. 

3. MONTECARLO SIMULATIONS 

In this section we use Montecarlo simulations to analyze 
the discrepancies observed by N09 and test the importance of 
radiation forces. We start from the observed distributions of 
continuum luminosity and broad line widths which, combined 
with assumed "real" distributions of bolometric correction and 
BLR column densities, allow us to obtain the real distribution 
of BH mass values and Eddington ratios. We consider two 
cases in which Mbh is obtained with or without important ra- 
diation pressure effect. By taking into account the intrinsic 
scatter in the Mbh ~ relation, we then estimate the observed 
distribution of stellar velocity dispersion from Mbh- Then we 
repeat the analysis by N09 and compare the "observed" dis- 
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tributions of BH mass and Eddington ratios derived from the 
different scaling relations for Mbh- 

For simplicity, we assume that the distributions of observed 
quantities and scaling parameters are lognormal. In particular, 
we assume that the observed distributions of L5100 and Vhp 
values can be expressed as 

iog(L5ioo)=/:=/:()+S£! (8) 

log(yff/3) = V = Vo + Sv7 (9) 

where / and j are normally distributed random numbers with 
zero average and unitary standard deviation. £0. M) are 
therefore the averages of the distributions of log(L5ioo) and 
log(V//^) values, while S£, Ey are the standard deviations. 
After the selection of the / and j random numbers, we apply 
the sample selection criteria adopted by N09 (42.8 < £+44 < 
44.8 and V > 1). The AGN are characterized by true distribu- 
tions of bolometric corrections and BLR cloud column densi- 
ties: 

log(/75ioo) = 6=So + Se/z (10) 
log(A?H)=AA = AAo + SA^^ (11) 

with the same conventions as above. It is then possible to 
recover the true distribution of BH masses from the £, V, B 
and M distributions. When radiation forces are important: 

log(MBH/M0) = = log (10^' + 10^^) (12) 

with 

Xl=(/rad + S/M) + 2V + 7/: (13) 

M2 = g+C + {21>-N) (14) 

where Y.fU represents the variance of the /rad scaling factor 
and we have assumed, for simplicity, that the variance in g 
is entirely dominated by the variance in the column density. 
Conversely, in the case where radiation forces are not impor- 
tant, the true distribution of BH masses is given by 

7Wv = (/+S/m) + 2V + 7/: (15) 

The observed distribution of stellar velocity dispersions is ob- 
tained, in both cases, from the Mbh - cr relation: 

Mk{qx Mm) = a + PS + Hmw (16) 

where S = log(o'/200kms"') and w is a normally distributed 
random number like /. Thus, it is possible to derive the "ob- 
served" BH masses and bolometric luminosities as: 

Mc = Mr^i{0vM^ir)-^MW (17) 

7Wvir = /+2V + 7£ (18) 

A^rad = log (10^"«'+2V+7C + 10^+^) (19) 

£boi = So + £ (20) 

where A^rad or represent the virial estimators with or 
without correction for radiation pressure. These provide the 
observed distributions of L/Lndd and, clearly, these observa- 
tional quantities do not include the true distributions of scaling 
parameter values. 

In fig. [T] we show the results for the case of BLR motions 
affected by radiation forces. We use the following set of val- 
ues: 5000 random realizations with (£o,S£) = (-0.2,0.3), 
(Vo,Sv) = (0-65,0-3), (60, Sg) = (0-95,0-2), (AAo,Eaa) = 
(23.0,0-5), = 0-3. This set of parameters was chosen to re- 
produce the AA„ and yMvir distributions observed by N09 but 
the actual adopted values do not influence our general conclu- 
sions. Our simulated samples nicely reproduces the features 
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Fig. 1. — Observed and true dtstributions of BH masses (top panels) and 
Eddington ratios (bottom panels) in the case of BLR motions affected by ra- 
diation forces. The red line denotes BH masses derived from the M^n — cr 
relation (JVdcr), the black line denotes BH masses computed using classi- 
cal virial relations (A^vir) and the blue line denotes BH_ masses computed 
using viiial relations corrected for radiation pressure (A^iaj). The shaded 
histograms indicate the true distributions of BH masses and Eddington ra- 
tios. The sharp cutoff in the observed L/Z-Edd distributions is present even if 
radiation pressure forces are important. 

of the A^rad and L/L^a distributions observed by N09: the 
distribution of BH masses computed with the radiation pres- 
sure correction is narrower than that based on Mbh - o" and 
classical virial relations; the observed distributions of L/Ledd 
ratios based on Mbh - cr and classical virial relations are in 
nice agreement, while that based on radiation pressure cor- 
rected virial masses is sharply truncated at L/Ledd — 0.15 and 
more sharply peaked. The true Mbh and L/Ledd distributions 
are intermediate between those found with and without the 
radiation pressure correction. For this simulation we have 
adopted Sjv^ = 0.5 dex but if we increase that value the true 
Mbh and L/L^dd distributions will be broader and approach- 
ing those observed when using Mbh without the radiation 
pressure correction. Allowing for an intrinsic dispersion in 
the scaling parameters, especially in A^h, smears the cutoff in 
the observed L/Lndd distribution which then approaches the 
true one. The comparison between the observed (blue line) 
and true (shaded histogram) L/Lndd distributions shows that a 
fraction of sources have Eddington ratios larger than the criti- 
cal value at which a BLR with A^h = 10^"' cm"^ becomes grav- 
itationally unbound. This is possible only because, allowing 
for a distribution of A^h values, a fraction of the sources have 
A^H > 10^^ cm"^. Then, if the true distribution of Eddington 
ratios extends to large values, BLRs are gravitationally bound 
only if their average column density increases with L/Ledd 
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Fig. 2. — Observed and true distributions of BH masses (top panels) and 
Eddington ratios (bottom panels) in the case of BLR motions NOT affected 
by radiation forces. Notation as in fig.[T] 

(see also |Dong et al.|2009| l. 

Finally, in fig. |2| we show the results for the case of BLR 
motions NOT affected by radiation forces. We use the same 
set of parameters as before, except for Sy = 0.25 to have sim- 
ilar true distributions of Mbh and L/Ledd- We observe again 
the same features that prompted N09 to consider radiation 
forces unimportant, and in particular the sharp cutoff in the 
observed L/Lndd distribution. In this case the true distribu- 
tions of Mbh and L/Ledd values are well matched by the ob- 
served ones, because the intrinsic scatter of the scaling pa- 
rameters is small compared to the combined scatter of the ob- 
served line widths and continuum luminosities. 

4. SUMMARY 



The simulations presented in figs. [T] |2| indicate that, when 
taken at face value, the comparison of the observed Mbh 
and L/L-Edi distributions would always lead to the conclu- 
sions that radiation forces are not important in determining 
the motions of BLR clouds. However, we have shown that it 
is not possible to assess the importance of radiation forces 
on BLR cloud motions on the basis of the observed Mbh 
and L/^Edd distributions. The differences between the L/LEdd 
distribution obtained using the Mbh - cr relation and the one 
based on virial masses with radiation pressure correction can 
also be explained by neglecting the intrinsic dispersion in the 
adopted scaling parameters (e.g. /rad and g). In particular, the 
sharp cutoff at L/Lndd — 0.15 observed by N09 corresponds 
to the critical luminosity at which radiation forces balance 
gravitational attraction on BLR clouds with the adopted A^h 
value (lO^-'cm"-') but does not indicate that radiation forces 
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are negligible. A broad distribution of A^h values (e.g. log- 
normal with lO^^cm"^ average and 0.5 dex standard devia- 
tion) will remove such sharp cutoff, fully explaining the ob- 
served differences in L/L-em distributions. The sources with 
L/L^dd > 0.15 are then those in the high tail of the A^h distri- 
bution. Moreover, at these large L/L^dd^ BLR clouds can be 
gravitationally bound only if they have large column densities, 
e.g. if A^H increases with L/L-^dd- In this scenario, the sources 
with A^H — 10^"^ cm"^, for which radiation pressure is negligi- 
ble, constitute a minority of the whole population. In conclu- 
sion, it is not possible to distinguish between the two scenar- 
ios in which radiation forces are negligible in all sources or in 
just a minority of them. 

When estimating BH masses using virial mass estimators, 
one should then always consider the possibility of important 
radiation forces by using the currently calibrated correction 
which corresponds to A^h — 10^"^ cm"^ (M08), until it is possi- 
ble to assess the possible dependence of A^h on the observed 



source properties. 

Finally, the comparison of the Mbh and L/Z-Edd distribu- 
tions proposed by N09 can still constrain the overall scatter 
in the scaling parameters used in virial mass estimators. A 
detailed analysis of this issue is beyond the scope of this pa- 
per but, briefly, after obtaining the true distributions of Mbh 
and L/L^dd from the observed distribution of stellar velocity 
dispersions one can match it with the Mbh and L/L^dd dis- 
tributions based on virial masses, convolved with a suitable 
broadening function. Such a broadening function would then 
provide the combined variance of BLR physical properties in 
the sample of AGN under examination. 
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